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Abstract 



We investigate the feasibility of analyzing deep inelastic structure functions in 
Hamiltonian formalism by combining the light-front BJL limit of high energy 
amplitudes and the Fock space (multi-parton) description of hadrons. This 
study is motivated by some of the theoretical questions emerging from the 
ongoing nonperturbative/perturbative studies in light-front QCD and also 
by current problems in the interface of perturbative/nonperturbative QCD. 
In this preliminary study, we address the unpolarized structure function F^. 
Our starting point is the expression for the quark structure function as the 
Fourier transform of the expectation value of the good component of a bilocal 
current in the target state. By expanding the target state in a set of multi- 
particle Fock states, the structure function is expressed as the sum of squares 
of multi-parton wavefunctions integrated over independent longitudinal and 
transverse parton momenta. Utilizing the fact that the multi-particle Fock 
states are connected with each other by the light-front QCD Hamiltonian, 
we study questions of cancellation of collinear singularities, factorization of 
mass singularities, and the logarithmic scaling violations in the Hamiltonian 
picture. In this paper the essential features of the formalism are illustrated 
utilizing the calculation of the structure function of a dressed quark and the 
evolution of the valence part of the qq bound state structure function up to 
order g 2 . 



*On leave of absence from Saha Institute of Nuclear Physics, Calcutta, India 
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I. INTRODUCTION 



One of the earliest motivations for studying field theories in the light-front formula- 
tion came from investigations in current algebra |IJ. Also of significant importance were 
the observation of Bjorken scaling in deep inelastic scattering and the development of 
Feynman's parton model 0. At a former level, starting from the algebra of currents in 
light- front theory and utilizing Bjorken- Johnson-Low (BJL) limit of high energy scattering 
one could "derive exact scaling" and arrive at an expression for the scaling function as the 
Fourier transform of a bilocal current matrix element j3J. The structure and consequences 
of light-front current algebra played a significant role in proposing Quantum Chromody- 
namics (QCD) as the underlying theory of strong interactions ||. However, many questions 
regarding renormalization aspects remained unanswered. 

Following the discovery of asymptotic freedom in QCD, scaling and scaling violations 
were studied || utilizing operator product expansions near the light cone and the theory 
of composite operator renormalization. In this language the emphasis is on manifest gauge 
and Lorenz covariance and the intuitive parton picture remains hidden. On the other hand, 
immediately after the discovery of scaling and prior to the discovery of asymptotic freedom, 
the utility of the infinite momentum frame picture/light-front language in providing sim- 
ple, intuitive description of high energy phenomena was recognized and exploited by many 
authors [[F|-[T0[|. Thus the intuitive approach (physical gauge, light cone variables, ...) to 



study various issues in deep inelastic scattering have been vigorously pursued [fn|-|I5] along 
with the operator product expansion formalism. In this approach the covariant Feynman 
Diagram (with four momentum conservation at vertices and off-mass shell particles in in- 
termediate states) is still employed. In contrast is the approach |p||T7jl using old-fashioned 
perturbation theory (three momentum conservation at the vertices, on-mass shell particles 
in intermediate states), light front gauge, and light front variables. As has been stressed 
many times by Brodsky, Lepage and collaborators [ 18H , it is the latter approach which keeps 
close contact with the intuitive parton model ideas within field theory. 

Old-fashioned perturbation theory in the infinite momentum frame was utilized by Drell, 
Levy, and Yan [7J to study scaling of the deep inelastic structure function in field theory 
in the context of pseudoscalar Yukawa model. After the emergence of QCD as the under- 
lying theory of strong interactions, various aspects of structure function have been studied 
in the old-fashioned light-front perturbation theory formalism. For example, 1) questions 



concerning the evolution of distributions with emphasis on end-point (x — > 1) behavior [17 



(higher twist effects), 2) phenomenological aspects of higher twist effects fI5[] , 3) gluon dis- 



tribution function p0| , 4) photon distribution function f21 |, and 5) theoretically motivated 
parameterizations of parton distributions []22j] . 

The main advantages of the Brodsky-Lepage framework may be summarized as follows. 
The light-front perturbation theory approach is the one closest to the intuitive parton pic- 
ture, while maintaining the field theoretic aspects of the problem in tact. The framework 
makes use of the properties of multi-parton wavefunctions which span both the perturbative 
and non-perturbative sector of the theory. The calculation of certain observables sometimes 
simplifies drastically compared to the other approaches to the problem. The utility of light 
front dynamics for the interpretation and the calculation of polarized structure functions 
and higher twist structure functions in general has been stressed also by Jaffe, Ji and collab- 
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orators ||23|| . The disadvantages of the approach are mainly associated with the complexities 
of renormalization in the Hamiltonian framework. (For early work on renormalization of 
QED in the infinite momentum frame see, for example, Ref. p4]j.) 

At present, two major approaches aimed at the non-perturbative study of field theories 
in the light-front formalism |25[ are the Discretized Light Cone Quantization (DLCQ) |26| 
and Light-Front Tamm-Dancoff (LFTD) approximation p7j. Recently there has also been 
a proposal to study QCD in the light front Hamiltonian approach with particular emphasis 



on the aspects of renormalization p8[ . Two of the most interesting features of light front 
dynamics in the Hamiltonian language are light front power counting P8,ZD and the special 
aspects of longitudinal dynamics. 

Perturbative and non-perturbative aspects of longitudinal dynamics are of interest both 
from theoretical and phenomenological considerations. The problem associated with exactly 
zero longitudinal momentum (the so-called "zero-mode" problem) has attracted lot of atten- 
tion recently. The divergences associated with very small longitudinal momentum also play 
a major role in the approach of Ref. [28[] in the context of incorporating non-perturbative 
effects in the effective light-front Hamiltonian via perturbation theory framework. Study 
of deep inelastic structure function in the language of multi-parton wavefunctions offers an 
opportunity to study directly the question of cancellation/non-cancellation of small longi- 
tudinal momentum divergences in physical observables. 

At small values of the longitudinal momentum (or equivalently at very small values of 
x), the long (longitudinal) distance aspects of the theory are probed and one enters the 
overlapping region of perturbative and non-perturbative dynamics. The small x behavior of 
structure functions has attracted lot of attention recently in the light of new data from the 



HERA facility. A recent work of interest is that of Mueller |30| on small x behavior using 
multi-parton light-front wavefunctions. 

In this work we study various aspects of the deep inelastic structure function F2 utilizing 
light-front symmetries and multi-parton wavefunctions. As is well known, some aspects of 
the problem are easily probed in coordinate space and some aspects are easily probed in 
momentum space. A convenient starting point to discuss features of deep inelastic scattering 
in the Hamiltonian framework is provided by the expression for the structure function as the 
Fourier transform of a bilocal matrix element. Among the many derivations of this result 
that exist in the literature, is the one utilizing equal-x + light-front current algebra and the 
BJL limit |4j]. In this approach (which closely follows the original derivation of Bjorken in 
the infinite momentum frame) partons are not introduced from the beginning and "perfect" 
scaling emerges (after canonical manipulations) as a result of special properties of the equal 
x + commutators of currents. One essential feature of the commutators is the non-local 
behavior in the longitudinal direction which is a consequence of the symmetries on the 
light-front. Furthermore, the current commutators are found to have the same form in free 
and interacting theories. Thus BJL limit together with light-front current algebra allows us 
to have a starting point (which is very close to the physical answer) for the Hamiltonian 
analysis of the structure function.. 

To address the problem of scaling violations in the Hamiltonian picture and the asso- 
ciated renormalization, it is convenient to introduce the constituents of the hadron. As is 
well known, one can make the synthesis between the current algebra picture and the parton 
picture by expanding the state in terms of multi-parton (constituent) wavefunctions pi 
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The partons that appear in this formalism need not be collinear or massless but they are al- 
ways on mass shell since we are using the Hamiltonian dynamics. The F 2 structure function 
appears as the sum of squares of multi-parton wavefunctions integrated over independent 
longitudinal and transverse momenta. The behavior of multi-parton wavefunctions is de- 
termined by the light-front Hamiltonian. In this picture scaling violations are intimately 
associated with the renormalization aspects of multi-parton wavefunctions. 

Apart from the issues of renormalization, there is a question that confronts any approach 
that attempts the description of structure function in terms of constituent wavefunctions. 
Is it possible, even in principle, to describe the hadron structure function in terms of a finite 
number of constituent wave functions, since " the average number is infinite" [0,0 according 
to the standard folklore. What we have to worry about are the gluon and sea distribution 
functions and their physical interpretation. To address this question, we need to take into 
consideration, their definition in terms of the multi-parton wavefunctions, various physical 
constraints from the normalization of the hadron state, longitudinal momentum sum rule, 
etc. The approach using parton wavefunctions offers the possibility to address this question. 

It is of interest also to see how various issues like the cancellation of collinear singularities, 
factorization (separation of soft and hard physical processes), suppression of coherent (inter- 
ference) effects in the hard processes, scale evolution of parton densities, etc in perturbative 
QCD emerge in the present formalism. 

The plan of this paper is as follows. In Sec. II we discuss the general feature of the 
structure function F 2 . Structure function of a dressed quark is discussed in Sec. Ill to 
illustrate the calculational framework. In Section IV we discuss the structure function of a 
meson-like bound state. Logarithmic scaling violations in the bound state structure function 
are discussed in Section V. Discussion and summary is presented in Sec. VI. Finally the 
details of the normalization of the state and longitudinal momentum fraction sum rules are 
discussed in appendices A and B respectively. 



II. UNPOLARIZED STRUCTURE FUNCTION F 2 

We start from the expression for the F 2 structure function (more precisely, the quark 
momentum fraction distribution function): 

F 2 (x) 



x 

'1 w 



J drje-^V^rj). (2.1) 

The generalized form factor V 1 is defined by 

(P\J + (y\0)\P) = P + V\ V ), (2.2) 

where 

J + (y I o) = ^{y)i + m - m^iy)}. (2.3) 

The bilocality is only in the longitudinal direction, i.e., y + = 0, y 1 - = 0. Note that x = —-p+ 
and r] = jP + y~. The field operator ifi + (y) is given by []33 



4 



where 

£(y) = £x.i I [dp x ][b{p 1)Sl )e-^- +d\p 1 ,-s 1 )ei^y-l ( 2 .5) 

si J 

with [dpi] defined in appendix A. In this framework, the operator has a simple structure and 
all the complexities are buried in the state | P). The state | P) can be expanded in terms of 
multi-parton Fock space amplitudes which are related to each other through the relativistic 
light-front version of the Schroedinger equation 

P- I P) = p\ } I P), (2-6) 

where P~ is the light-front QCD Hamiltonian which can be written as 

P~ = P +V (2.7) 

with Pq is the free part and V is the interaction part. The explicit expressions for P~ and 
V are given in Ref. ||33|| . We note that the main focus of the ongoing effort in light-front 



theory is to solve this equation nonperturbatively. We can explore various features of the 
structure function utilizing the light-front Hamiltonian P~ . 

When we address the problem of calculating the structure function of a composite system 
like a baryon or a meson, which of course are the physically interesting systems, we are 
immediately bombarded with all the complexities in the real world. To break up the problem 
in to simpler pieces, in the following section we first consider an artificial target, namely, 
a dressed parton. Even though such an object does not appear in the real world, this 
calculation will help us to understand some of the key issues in a simpler setting and also to 
set the notation. For calculations of structure function of a quark done in the framework of 
Feynman perturbation theory in planar gauge see Dokshitzer, Dyakonov, and Troyan [13]. 



III. STRUCTURE FUNCTION OF A DRESSED QUARK 

A. Parton picture 

We take the state | P) to be a dressed quark and expand this state in terms of bare 
states of quark, quark plus gluon, quark plus two gluons, etc. The expansion takes the form 

|P<r) = &(i»&t(P,<7)|0) 

+ E / {*i} / W 

ci,A2 



2 (P, a | h, m; k 2 , \ 2 )p{2nfP+5 z {P - h - k 2 )b\k u a 1 )a\k 2 , A 2 ) | 0) 
+ \ E J{dh} J{[dk 3 } 



o"i,A2,A3 ' 

f3/ 



fc(P, a | hi, <7i; k 2 , A 2 ; k 3 , \ 3 )^2(2n)zp+5 s (P ~h-k 2 - k 3 ) 

b\k u al)a\k 2j \ 2 )a\k 3 ,X 3 ) | 0) 
+ .... (3.1) 
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The factor of ~ in front of the expression in the third term above is the symmetry factor 
for identical bosons. 

The function <ft(P, a) is the probability amplitude to find a bare quark with momentum 
P and helicity o in a dressed quark, the function <f)(P,<J \ &iO"i, &2A2) is the probability 
amplitude to find a bare quark with momentum k\ and helicity e\ and a bare gluon with 
momentum /c 2 and helicity A 2 in the dressed quark, etc. 

We evaluate the expression in Eq.( |2.1| ) explicitly, noting that in the present case the 
contribution from the second term in Eq. (|2.3j ) is zero. Introduce the Jacobi momenta (xi, n^) 

— XiP~^~ , k^ = -\- XiP (3.2) 

so that 

]>> = 1,$X = 0. (3.3) 



Also introduce the amplitudes, 



01 = fa, 

UK,ki) = ^M^4), (3-4) 



and so on. We get 

F 2 (x) 



x 



+ ^2 J dx 2 J d 2 K^ J d 2 K,25(l — x — X2)5 2 {k^ + k 2 ) I ip^'^ix, x 2 , n 2 ) 
+ - J dx 2 J dx 3 J d 2 K^ J d 2 K 2 J d 2 K%5(l — x — x 2 — x^)5 2 {k{ 



-L 1 „-L-\ 
2 



fl,A2,A3 ' 



^ lA2A3 (x,^;x 2 ,4;x 3 ,4) I 2 + ... (3.5) 



This equation makes manifest the partonic interpretation of the quark distribution function, 
namely, the quark distribution function of a dressed quark is the incoherent sum of prob- 
abilities to find a bare parton (quark) with longitudinal momentum fraction x in various 
multi-particle Fock states of the dressed quark. Since we have computed the distribution 
function in field theory, there are also significant differences from the traditional parton 
model . Most important difference is the fact that the partons in field theory have trans- 
verse momenta ranging from zero to infinity. Whether the structure function scales or not 
now depends on the ultraviolet behavior of the multi-parton wavef unctions. By analyzing 
various interactions, one easily finds that in superrenormalizable interactions, the transverse 
momentum integrals converge in the ultraviolet and the structure function scales, whereas 
in renormalizable interactions, the transverse momentum integrals diverge in the ultraviolet 
which in turn leads to scaling violations in the structure function. 

We should also note that even though (|3.5|) looks like an incoherent sum, interference 
effects are also present in this expression. We will elaborate more on this in Sec.V. 
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B. Dressing with one gluon 



In this section we evaluate the structure function of a dressed quark to order a s . Our 
starting point is the eigenvalue equation obeyed by the dressed quark state, i.e., Eq. Q2.6Q. 
We substitute the Fock space expansion ( |3.1|) in (|2.6| ) and make projection with a bare one 
quark - one gluon state tf(p\, si)a^(j>2, P2) | 0) keeping only terms up to order g. We arrive 
at 

M 2 + (P x ) 2 m 2 + (pi) 2 (p^) 2 



p+ pi pi 



=T a ^=xU 2 ZT t^-tr 1 — ]x ff .(e£) (3-6) 



2(2tt) 3 y/p+ Pi P + 

We rewrite the above equation in terms of Jacobi momenta (pi = XiP + , &i = pi + XiP^) 
and the wavefunctions fa which are functions of Jacobi momenta. Using the notation x = 
Xi, K\ — k and using the facts X\ + x 2 = 1, K\ + k 2 = 0, we have 

r7ix2 m 2 + (k^) 2 (k l ) 2 , B1P 21, _|_ 1 J_x 

M fore ;1-x,-k ) = 

X 1 — X 

5 ™ 1 t r o K± o^-.k 1 - — im ± ± . ± * 

~ l^f 7r=^ xU " 2 ~ i — a ~ a im]xA ^ ] * (3 ' 7) 

Taking the bare and dressed quarks to be massless, we arrive at 

E / ^ I ^fore^l-*,-^) \ 2 = -£- 3 C f I fa I 2 / d2 ^Ay ( 3 - 8 ) 



where C/ = 2 n ^ e t ransverse momentum integral is divergent at both limits of integra- 
tion. We regulate the lower limit by \i and the upper limit by Q. Thus we have 

— =| fa I [5(1 - x) + - Cf —l n -\. (3.9) 
The normalization condition reads 

1 ' 2 ...Q 

Within the present approximation (valid only upto a s ) 



fa | 2 [1 + ^C f [ dx\±^-ln%] = 1. (3.10) 



a s r l + x 2 Q 2 

— Of / ax In— 

2tt 'J l-x fi 2 



fa | 2 = 1 - -±C f / dx— In^r. (3.11) 



In the second term we recognize the familiar expression of wavefunction correction of the 
state n in old fashioned perturbation theory, namely, Y,L Ir^rl ■ 

III' yJl/ n &m) 

Thus to order a s , 

— = 5(1 - x) + _C>1„_[_ - 5(1 - x) y dy— ]. (3.12) 
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Note that (|3.12|) can also be written as 



_ = *(1 - *) + ^^^[jjT^ + 2 5(1 " ^ (3 ' 13) 

which is a more familiar expression. Note that by construction | -02 ( probability 
density. However, this function is singular as x — > 1 (gluon longitudinal momentum fraction 
approaching zero). To get a finite probability density we have to introduce a cutoff e (x g i uon > 
e), for example. In a physical cross section, this e cannot appear and here we have an 
explicit example of this cancellation. The probabilistic interpretation of the splitting function 
Pqq = ~[37 which arise from real gluon emission is obvious in our derivation. On the other 
hand, the function P qq = + 1^(1 — x ) doesn't have the probabilistic interpretation 

since it includes contribution from virtual gluon emission. This is immediately transparent 
from the relation 

" dxP qq (x) = 0. (3.14) 

We also note that the divergence arising from small transverse momentum (the familiar mass 
singularity) cannot be handled properly in the present calculation. This is to be contrasted 
with the calculation of the meson structure function (see the following section) where the 
mass singularities can be properly absorbed into the nonperturbative part of the structure 
function. 

It is instructive to compare the above derivation with typical calculations that exist 
in the literature. Conventionally, one evaluates the imaginary part of the forward virtual 
Compton amplitude order by order in perturbation theory. In the lowest order, in addition 
to the real gluon emission, one also has various self-energy diagrams and the dressing of 
the photon vertex by the gluon. The ultraviolet part of the vertex and parts of the self- 
energy contributions cancel each other as a result of the QED Ward Identity Z\ = Z 2 . 
Thus, for example, in planar gauge, Dokshitzer, Dyakonov and Troyan |L3| find by explicit 



calculations that neither collinear nor infrared mass singularities affect the "partonometric 
vertex" and that the radiative corrections to the photon vertex appear to be effectively unity. 
By explicit calculations both in the 3 theory and QCD, Collins, Soper, and Sterman p£3 



find that non-ladder diagrams are either higher twist or contribute only to the hard part. 
In the calculation of Brodsky and Lepage [17] using light-front time ordered perturbation 



theory in the gauge A + = 0, the identity Z\ = Z 2 is used. In the present calculation, since 
we have started from the BJL limit, dressing of the photon vertex, gluon radiation from the 
final state, and parts of self-energy contributions (that cancel eventually) does not appear 
at all which leads to simplification of the calculation. It is of interest to see whether such 
simplifications persist to higher orders in the present framework. 

IV. STRUCTURE FUNCTION OF A MESON-LIKE BOUND STATE 

For simplicity we consider a meson-like bound state. We expand the state | P) for qq 
bound state in terms of the Fock components qq, qqg, ... as follows. 



S 



P) = E / [Mi] J [dh 



fl,CT2 



2 (P I kx,ax; k 2 , a 2 )^2((27r) 3 P+5 3 (P - fej - k 2 p(h, a 1 )S(k 2 , a 2 ) | 0) 
+ E / {*i> / {*2> / {dk 3 } 



C1,(T2,A3 ' 



3 (F | fci, a i; fc 2 , a 2 ; fc 3 , \ 3 )^2(2tt) 3 P+8 3 (P - A* - k 2 - k 3 ) 

b\k 1 ,a 1 )d\k 2 ,a 2 )a\k 3 ,X 3 ) | 0) 
+ ... . (4.1) 

Here <fr 2 is the probability amplitude to find a quark and an antiquark in the meson, is 
the probability amplitude to find a quark, antiquark and a gluon in the meson etc. 

As in Sec. Ill we evaluate the expression in Eq. (|2.1|) explicitly. The contribution from 
the first term (from the quark), in terms of the amplitudes 

MK,kt) = ^:Mxi,Ki), (4-2) 



and so on, is 

F 2 (x)« 



X 



J dx 2 J d 2 n^ J d 2 K 2 8(l — x — x 2 )5 2 (ki + k 2 ) \ il) 2 1,cr2 {x, k^] x 2 k 2~) | 2 

0"1,(T2 

+ J dx 2 J dx3 J d 2 K^ J d 2 K 2 J d 2 K^5(l — x — x 2 — xz)5 2 {ki + k 2 + K3) 



0-l,0"2, A3 ' 



r 3 ua2M (x,Ki;x 2 ,4;x 3 ,Ki) | 2 +... . (4.3) 



Again, the partonic interpretation of the F 2 structure function is manifest in this expression. 
Contributions to the structure function from the second term in Eq.( [2.1| ) is 

2 ^ - = E / d x 2 [ d 2 K^ I d 2 K 2 5(l — x — x 2 )S 2 (ki + k 2 ) \ i/j 2 1,cr2 (x 2 , k 2 ; x, k^) \ 2 

+ ^2 J dx 2 J dxs J d 2 &i J d 2 n 2 J d 2 n^6(l — x — x 2 — x 3 )5 2 {ni + k 2 + k 3 ) 

<71,<T2>A3 

I ^ 2M {x 2l 4; x, 4; x 3 , 4) | 2 +... . (4.4) 
The normalization condition given in Eq.( [B3"D guarantees that 

J x x 

which reflects the fact that there are two valence particles in the meson. Since the bilo- 
cal current component J + involves only fermions explicitly, we appear to have missed the 
contributions from the gluon constituents altogether. Gluonic contribution to the struc- 
ture function F 2 is most easily calculated by studying the hadron expectation value of the 
conserved longitudinal momentum operator P + . The details are given in appendix C. 
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From the normalization condition, it is clear that the valence distribution receives contri- 
bution from the amplitudes ^ "03, ••• at any scale /jl. This has interesting phenomenological 
implications. In the model for the meson with only a quark-antiquark pair of equal mass, the 
valence distribution function will peak at x = ~. If there are more than just the two particles 
in the system, the resulting valence distribution will naturally have an enhancement in the 
x < g region and a depletion in the x > | region, as a simple consequence of longitudinal 
momentum conservation. 

The equation (|4.3|) as it stands is useful only when the bound state solution in QCD 
is known in terms of the multi-parton wavef unctions. The wavefunctions, as they stand, 
span both the perturbative and non-perturbative sectors of the theory. Great progress in 
the understanding of QCD in the high energy sector is made in the past by separating 
the soft (non-perturbative) and hard (perturbative) regions of QCD via the machinery of 
factorization. It is of interest to see under what circumstances a factorization occurs in the 
formal result of Eq. (|4.3| ) and a perturbative picture of scaling violations emerges finally. 
We address this issue in the following section. 



V. PERTURBATIVE PICTURE OF SCALING VIOLATIONS 

To address the issue of scaling violations in the structure function of the "meson-like" 
bound state, it is convenient to separate the momentum space into low-energy and high- 
energy sectors. Such a separation has been introduced in the past in the study of renor- 
malization aspects of bound state equations |Q in light-front field theory. The two sectors 



are formally defined by introducing cutoff factors in the momentum space integrals. How 
to cutoff the momentum integrals in a sensible and convenient way in light-front theory 
is a subject under active research at the present time. Complications arise because of the 
possibility of large energy divergences from both small k + and large k regions. In the fol- 
lowing we investigate only the effects of logarithmic divergences arising from large transverse 
momenta, ignore subtelities arising from both small x (x — > 0) and large x {x — > 1) regions 
and subsequently use simple transverse momentum cutoff. For complications arising from 
x — > 1 region see Ref . [TL7] . 



A. Scale separation 

We define the soft region to be k < fi and the hard region to be [i < k < A. /i serves 
as a factorization scale which separates soft and hard regions. Since it is an intermediate 
scale introduced artificially purely for convenience, physical structure function should be 
independent of /i. The multi-parton amplitude ip% is a function of a single relative transverse 
momentum nj; and we define 

The amplitude ip 3 is a function of two relative momenta, k± and and we define 
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V>3 



^| s , < 4,4 < ,U 
l/jf, < 4 < /I, 
^3 S , ii < 4 < A, < 4 < 11, 

i>% h , fi < 4,4 < a 



(5.2) 



Let us consider the quark distribution function q{x) = F2 ^ x ' defined in Eq. fl4.3p . In presence 
of the ultraviolet cutoff A, q(x) depends on A and schematically we have, 



g(x,A 2 ) = ErV 2 2 + E/ A 

Jo Jo Jo 



For convenience, we write, 



q(x,A 2 ) = q 2 (x,A 2 )+q 3 (x,A 2 ). 



(5.3) 



(5.4) 



where the subscripts 2 and 3 denotes the two-particle and three-particle contributions re- 
spectively. Schematically we have, 



q(x,A 2 ) = q(x,n 2 ) + Y J / IV> 



h |2 
2 I 



1 2 



Jfl 

A /-A 



A rp 



H JO 



1>. 



hs 1 2 



W | 2 



(5.5) 



We investigate the contributions from the amplitudes and ?/>3 S to order a s in the follow- 
ing. 



B. Dressing with one gluon 

We substitute the Fock expansion Eq. ( [Lip in Eq. (2~TS) and make projection with a 
three particle state 6^(A;i, 0"i)g4&2j er 2 )a^(A;3, cr 3 ) | 0) from the left. In terms of the amplitudes 
^2, ^3, we get, 

^ a2X3 (x, «i; x 2 , k 2 ; 1 - x - x 2 , k 3 ) = -Mi + -M 2 , (5.6) 

where 

Mi = r= r " /1 1 ^ ^(1 -x 2 ,-4;x 2 ,4), (5.7) 

^ y/2(2^) 3 V1-X-X2 



and 



M 2 = i^=L==T a ^==L— y 2 ^(*, «r; 1 - x, -4) (5.8) 



with 
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E = [M 2 — 



, m 2 + (k^) 2 m 2 + (k^) 2 
x x 2 



^2 



.x 



and 



14 = ^11,-^ g i (-^-M H „. 

^ 1 - X - X 2 

°2 



1 — X — £ 2 

(0-^.^2 — im) 
1 - x 2 

cr 1 ./^ — im) 



(5.9) 
(5.10) 



y 



l-x 



C. Perturbative analysis 

For nj- hard and soft, ftf- + ~ «i" and the multiple transverse momentum integral 
over -03 factorises into two independent integrals and the longitudinal momentum fraction 
integrals become convolutions. The contribution from M. x to ^3 is, 



^3,l' Cr2 ' A3 ( X > X 2i K 2~i 1 — X — X 2 , — K 2~) — / „ = ; = ^ 



™xVl-x-a;2 1 



^2(2^)3 1 " ^2 (^) 2 

2k£ h / / j 



Thus the contribution from .Mi to the structure function is 

,2 



where 



^'^(l -£2,-^; £2,/^). ( 5 - 12 ) 
function is 

E/l^r-^^jffw^A (5.i3) 

x l + (-) 2 

^) = T 3^- (5-14) 

y y 

For the configuration k± hard, soft, contribution from M. 2 does not factorize and the 
asymptotic behavior of the integrand critically depends on the asymptotic behavior of the 
two-particle wavefunction -0 2 . To determine this behavior, we have to analyze the bound 
state equation which shows that for large transverse momentum r ip2{^' L ) ~ ^xp ■ Thus 
contribution from M 2 for scale evolution is suppressed by the bound state wavefunction. 
Analysis of the interference terms (between Mi and M 2 ) shows that their contribution also 
is suppressed by the bound state wavefunction. 

For the configuration soft, hard, contributions from A4i and the interference terms 
are suppressed by the wave function. Contribution from Ai 2 factorises both in transverse 
and longitudinal space and generate a pure wavefunction renormalization contribution: 

E/i«r=|^/^, 2 (,,^). (5.i5) 

We have seen that even though the multi-parton contributions to the structure function 
involve both coherent and incoherent phenomena, in the hard region coherent effects are 
suppressed by the wavefunction. 
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D. Corrections from normalization condition 

In the dressed quark calculation, we have seen that the singularity that arises as x — > 1 
from real gluon emission is canceled by the correction from the normalization of the state 
(virtual gluon emission contribution from wave function renormalization). In the meson 
bound state calculation, so far we have studied the effects of a hard real gluon emission. In 
this section we study the corrections arising from the normalization condition of the quark 
distribution in the composite bound state. 

Collecting all the terms arising from the hard gluon emission contributing to the quark 
distribution function, we have, 

q(x, A 2 ) = q 2 (x, fi 2 ) + q 3 (x, fi 2 ) 

a Sn . A 2 f 1 dy x 2 

A 2 

+ ^C f ln-q 2 (x,^ 2 ) J dyP(y). (5.16) 

We have a similar expression for the antiquark distribution function. 

The normalization condition on the quark distribution function should be such that there 
is one valence quark in the bound state at any scale Q. We choose the factorization scale 
/jl — Q . Let us first set the scale A = Q . Then we have (in the truncated Fock space) 

f 1 dxq 2 (x, Ql) + f 1 dxq 3 (x, Q 2 ) = 1. (5.17) 

JO JO 

Next set the scale A = Q. We still require 

C dxq 2 (x, Q 2 ) + C dxq 3 (x, Q 2 ) = 1. (5.18) 

JO JO 

We note that the evolution of q 3 requires an extra hard gluon which is not available in the 
truncated Fock space. Thus in the present approximation q 3 (x, Q 2 ) = q 3 (x, Ql). 
Carrying out the integration explicitly, we arrive at 

fi 2a O 2 r r 1 

j o dxq 2 (x,^ 2 )[l + -^-C f ln^ J dyP(y)} + jf dxq 3 (x,Q 2 ) = 1 (5.19) 

Thus we face the necessity to "renormalize" our quark distribution function. Let us 
define a renormalized quark distribution function 

tf(x, Ql) = q 2 (x, Ql) [1 + ^C f ln^ / dyP(y)\ (5.20) 

Z7T Qq JO 

so that, to order a s , 

f 1 dxqX(x, Ql) + / dxq 3 (x, Q 2 ) = 1. (5.21) 

JO JO 

We have, 
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n Cr r 1 

q 2 (x,Cft) = q*(x,<%)[l-2^Cfln% i dyP(y)\. (5.22) 

Collecting all the terms, to order a s , we have the normalized quark distribution function, 

q(x,Q 2 )=q?(x,Ql) 

a O 2 r 1 r 1 
+ 7T C f ln 7v dyq«(y,Q 2 ) dz5(zy - x)P(z) 
2,ti Jo Jo 

+ q 3 (x,Q 2 ) (5.23) 

with P(z) = P{z) - 5(z - 1) Jo 1 dyP(y). 

We see that just as in the dressed quark case, the singularity arising as x — > 1 from 
real gluon emission is canceled in the quark distribution function once the normalization 
condition is properly taken in to account. From this derivation we begin to recognize the 
lowest order term of the Altarelli-Parisi evolution equation. 



VI. DISCUSSION 

In this paper we have investigated the feasibility of analyzing deep inelastic structure 
functions in the Hamiltonian formulation. Our starting point has been the non-perturbative 
result obtained using equal-x + light-front current algebra and the BJL limit which clearly 
shows exact Bjorken scaling if we forget the non-trivial renormalization issues involved in the 
matrix element of the bilocal current sandwiched between hadronic states. In this prelimi- 
nary study we have focused on the leading contribution (ignoring power suppressed terms) 
to the unpolarized structure function F2. Since the good (+) component of the bilocal cur- 
rent is involved, expanding the hadronic state into multi-parton wavef unctions, we obtained 
the structure function — as a sum of the squared multi-parton wavefunction integrated 
over longitudinal and transverse momenta. Taking into account the QCD interaction which 
connects the various wavefunctions and scaling violation comes in the picture. We took the 
light-front Hamiltonian and demonstrated how one could get logarithmic growth and the 
splitting functions in the context of a dressed quark. The probabilistic interpretation of the 
splitting function is very clear in this language. We also noted how x — > singularity could 
be taken care of by the normalization of the state in this language. In the bound state cal- 
culations involving a meson, we clearly illustrated how factorization of the non-perturbative 
and perturbative parts in the structure function occurred in our approach. We have also 
illustrated how various coherent effects are power suppressed in the hard region which leads 
to the standard perturbative evolution of parton distributions. 

It is of interest to carry out higher order calculations for the dressed quark to see explicitly 
how the running coupling constant comes into play in this formulation. 

The ultimate test of the present formulation (success or failure) is in its application to 
the analysis of the so-called higher-twist structure functions. The longitudinal structure 
function Fl and the polarized structure function Gt currently attracts lot of attention. 
Even the questions of physical interpretation, factorization and evolution of these structure 
functions are under intense investigation and debate at the present time. We plan to study 
these problems in the near future in the formulation presented in this paper which keeps 
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close contact with physical intuitions of the pre-QCD parton picture while maintaining the 
complexities of QCD. 
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APPENDIX A: NOTATIONS AND CONVENTIONS 

We use the definitions 

x ± = x°± x 3 , x x = (x\ x 2 ). (Al) 
7 ± = 7°± 7 3 . (A2) 

^± = A ± -0, withA ± = ^7 T 7 ± - (A3) 
The normalization of the state is 

(P | P') = 2{2irfP + 5\P - P'). (A4) 

The volume elements 

= dk + <Pk+ 
1 J 2(2^)3^ 

<-,-.-. d k , . . 

{dk} = , (A6) 



^2(2^)^+ 

and 
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APPENDIX B: NORMALIZATION OF STATE 

The normalization of the state | P) is 

(P' | P) = 2(2n) 3 P + S 3 (P - P'). (Bl) 
In the truncated space for the meson state, we have, 

]T f dk+d 2 k^ J dk+d 2 k 2 L 5 3 {P -h- k 2 ) | (j) 2 {P | k 1} (7i; k 2 , a 2 ) | 2 
+ E / dktd*k£ J dk+d 2 k 2 L dk+d 2 k^5 3 {P -k ± -k 2 - h) 

x | 3 (P | k u (7i ; fc 2 , a 2 , k 3 , A3) | 2 = 1. (B2) 



01,(72 ' 



C1,(T2,A3 ' 



In terms of the amplitudes ip2,i/>3, the normalization condition in the truncated Fock space 
sector reads 

2 y dzidV y da; 2 rf 2 45(l - xi - :r 2 )5 2 (4 + 4) | ^'^i, x 2 , 4) | 2 + 
y dx\d 2 Ki J dx 2 d 2 K2 J dxsd 2 K^S(l — x\ — x 2 — x 3 )5 2 (4 + 4 + 4) 

x 1 ^ 1,CT2,A3 (zi,4;z2,4;z3,4) l 2 = i- (B3) 



C1,<T2 



APPENDIX C: OPERATOR DEFINITION OF GLUON DISTRIBUTION 
FUNCTION AND THE LONGITUDINAL MOMENTUM SUM RULE 

Consider the conserved longitudinal momentum operator 

P+ = I y dx-d 2 x ± 9 ++ (CI) 

where 

0++ = 0++ + 0++. (C2) 

The fermion contribution 

#++ = ifa+d+ij) (C3) 

and the gluon contribution 

6>++ = p +va F+ a (C4) 

where 

_ Q^A ua — d u A^ a + g f abc A vc (C5) 

From the definition 
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Ff(x) = ^ J d V e-^(P | myh + m - miMv)] I P), (C6) 
with r\ = ^P + y~, we have 

/ dxFf(z) = (-^) 2 <P | e + F + I P). (C7) 
Formally, we can define the "gluon structure function" |34| 

F 2 G (x) = ^ J dy-e-^-^P | F + » a (y-)F+ a (0) \ P). (C8) 

We have, 

JdxF 2 G (x) = ( 1 L) 2 (P\e+ + \P} (C9) 
and the momentum sum rule 

J dx[F 2 F + if] = 1. (CIO) 

Next we explicitly verify the longitudinal momentum sum rule in the truncated Fock 
space. The fermionic part of the longitudinal momentum operator 

P+ = - J dx~d 2 x ± 2i(^j + yd + ^ + . (Cll) 

Using the field expansion 

^ + (*) = E Xx [ [dp][b(p, A)e"^ + <fl(p, -\)e^% (C12) 

A J 

K = E / (dp)p + [bHp, A)6Cp, A) + d\p, A)d(p, A)], (C13) 

with (dp) defined in appendix A. 

In A +a = gauge, the gluonic part of the longitudinal momentum operator 

p+ = - j dx~d 2 x ± d + A ja d + A ja . (C14) 
2 J 

Using the field expansion 

A ja = E / (dq)[e{a a (q, A)e"^ + (e^)* K)t(g, \)e^% (C15) 
a J 

P g = E / (^)<? + K)% A)a a (g, A). (C16) 
a J 

We have the relation 
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(P' I P + I P) = [2(2tt) 3 P + 5 3 (P - P')]P + . (C17) 
Explicit evaluation in the truncated Fock space gives 

J dxid 2 K^ J dx 2 d 2 ^5(1 — X\ — x 2 )5 2 (k^ + 



K 2 

x (xi + x 2 ) I ip% 1,<T2 (xi, re^; x 2 , «2 ) I 2 + 



J 1 43 



y dx\d 2 Ki J dx 2 d 2 K 2 J dx 3 d 2 K^ 5(1 — x\ — x 2 — x^& 2 (k^ + k 2 + 

x (xi +x 2 + x 3 ) | ^3 1 ' CT2 ' A3 (xi, 4; x 2 , x 3 , nj) | 2 = 1, (C18) 



Cl,^2,A3 



which is automatically satisfied because of the normalization condition given in Eq . (|B3| ) . 
Thus momentum fraction sum rule is trivially satisfied given the normalization condition on 
the amplitudes ifj 2 ,ip3 etc. 



18 



REFERENCES 



[1] S. ADLER and R. DASHEN, "Current Algebra and Applications to Particle Physics", 
Benjamin, New York, 1968; V.DE ALFARO, S. FUBINI, G. FURLAN, and C. ROS- 
SETTI, "Currents in Hadron Physics", North- Holland Publishing Company, Amster- 
dam, 1973; R. JACKIW, in "Current Algebra and Anomalies", (S. B. Treiman et al, 
Eds.) Princeton University Press, Princeton, New Jersey, 1985. 

[2] J.D. BJORKEN, Phys. Rev. 179 (1969) ,1547. 

[3] R.P FEYNMAN, "Photon-Hadron Interactions", Benjamin, Reading, MA, 1972. 
[4] R. JACKIW, in Springer Tracts in Modern Physics, 62, (1972). 

[5] H. FRITZSCH, M. GELL-MANN, and H. LEUTWYLER, Phys. Lett. B 47 (1973), 365 
and references therein. 

[6] H. GEORGI and H.D. POLITZER, Phys. Rev. D 9 (1973), 416 ; D.J. GROSS and F. 

WILCZEK, Phys. Rev. D 9 (1973), 980. 
[7] S.D. DRELL, D.J. LEVY, and T.M. YAN, Phys. Rev. 181 (1969). 2159; Phys. Rev. D 

1 (1970), 1617. 

[8] J.D. BJORKEN, J.B. KOGUT, and D.E. SOPER, Phys. Rev. D 3 (1971), 1382. 

[9] J.B. KOGUT and L. SUSSKIND, Phys. Rep. C 8 (1973), 75. 
[10] L.N. LIPATOV, Yad. Fiz. 20 (1974), 181 [Sov. J. Nucl. Phys. 20 (1975), 94]. 
[11] D.E. SOPER, Phys. Rev. D 15 (1977), 1141. 

[12] G. CURCI, W. FURMANSKI and R. PETRONZIO, Nucl. Phys. B 175 (1980), 27. 
[13] YU.L. DOKSHITZER, D.I. DYAKONOV, and S.I. TROYAN, Phys. Rep. C 58 (1980), 
269. 

[14] J.C. COLLINS and D.E. SOPER, Nucl. Phys. B 194 (1982), 445. 

[15] A.H. MUELLER, Phys. Rep. C 73 (1981), 237. 

[16] K.J. KIM and K. SCHILCHER, Phys. Rev. D 17 (1978), 2800. 

[17] S.J. BRODSKY and CP LEPAGE in "Quantum Chromodynamics" , (A. Mosher, Ed.), 
Proceedings of Summer Institute on Particle Physics, July 9-20, 1979, SLAC Report No. 
224, (1980). 

[18] S.J. BRODSKY and CP LEPAGE, in "Perturbative Quantum Chromodynamics", (A. 

Mueller, Ed.), World Scientific, Singapore, 1989, and references therein. 
[19] J.F. GUNION, P. NASON, and R. BLANKENBECLER, Phys. Rev. D 29 (1984), 2491 

and references therein. 

[20] S.J. BRODSKY and J.F. GUNION, Phys. Rev. D 19 (1979), 1005; S.J. BRODSKY and 

I. SCHMIDT, Phys. Lett. B 234 (1990), 144. 
[21] M. BURKARDT, Nucl. Phys. B 373 (1992), 371. 

[22] S.J. BRODSKY, M. BURKARDT I. SCHMIDT, Nucl. Phys. B 441 (1995), 197. 

[23] R.L. JAFFE, "Spin, Twist, and Hadron Structure in Deep Inelastic Scattering", MIT 

preprint MIT-CTP-2506, January 1996, and references therein. 
[24] S.J. BRODSKY, R. ROSKIES and R. SUAYA, Phys. Rev. D 8 (1973), 4574. 
[25] See for example, "Theory of Hadrons and Light-Front QCD", (St. D. Glazek, Ed.), 

World Scientific, Singapore, 1995 and references therein. 
[26] H.C. PAULI and S.J. BRODSKY, Phys. Rev. D 32 (1985), 1993; ibid. (1985), 2001. 
[27] R.J. PERRY, A. HARINDRANATH, and K.C WILSON, Phys. Rev. Lett. 65 (1990), 

2959; R.J. Perry and A. Harindranath, Phys. Rev. D 43 (1991), 4051. 



19 



[28] K.G. WILSON, T.S. WALHOUT, A. HARINDRANATH, W.-M. ZHANG, R.J. 

PERRY, and St. D. GLAZEK, Phys. Rev. D 49 (1994), 6720. 
[29] A. HARINDRANATH, "Consequences of scale invariance and scale symmetry breaking 

in light-front dynamics", preprint, March 1996. 
[30] A. H. MUELLER, Nucl. Phys. B 415 (1994), 373. 

[31] H. LEUTWYLER, in "Broken Scale Invariance and the Light Cone", (M. Dal Cin, G.J. 
Iverson, and A. Perlmutter, Eds.), Gordon and Breach Science Publishers, New York, 
1971. 

[32] J.D. BJORKEN in "QCD Hard Hadronic Processes", (Bradley Cox, Ed.), Plenum Press, 
New York, 1988. 

[33] W-.M. ZHANG and A. HARINDRANATH, Phys. Rev. D 48 (1993), 4903. 

[34] J.C. COLLINS, D.E. SOPER, and G. STERMAN, in "Perturbative Quantum Chro- 

modynamics", (A. H. Mueller, Ed.) World Scientific, Singapore, 1989 and references 

therein. 

[35] St. D. GLAZEK, A. HARINDRANATH, S. PINSKY, J. SHIGEMITSU, and K.G. WIL- 
SON, Phys. Rev. D 47 (1993), 1599. 



20 



